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Telephone structures must remain in operation during the shock and 
vibration caused by blasts, earthquakes, or other dynamic forces. Tradi- 
tionally, various numerical methods, including the finite element approach, 
have been used to analyze such problems. These methods, although effective, 
generally require excessive and costly computations. In contrast, the Fourier 
transform method used in conjunction with a fast Fourier transform 
algorithm is much more economical. In addition, shock and vibration 
problems involving frequency-dependent parameters can be effectively 
treated by the Fourier method. However, to make the Fourier method more 
effective and widely applicable, various tractable input-transfer-output 
relations are needed. This paper derives a set of simple equality and in- 
equality relations that allow various response parameters to be con- 
veniently estimated based on partial knowledge of the input and the 
structures. In particular, two lower bounds and eight upper bounds of the 
maximum response of linear structures are presented. Simple structures 
subjected to impulse loads resembling blasts and a random transient load 
resembling earthquakes are studied. Practical applications of the method 
to telephone structures are demonstrated by the analysis of a battery stand 
and a community dial office system in an earthquake area. 

I. INTRODUCTION 

The Fourier transform method has long been used to characterize 
linear systems and to identify the frequency content of waveforms. 
Sneddon 1 demonstrates classical applications of the Fourier transform 
method to engineering and physics problems. Amba-Rao 2 gives recent 
applications of this method to elasticity, White 3 illustrates system 
response calculations, Le Bail 4 demonstrates boundary-value problems 
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in physics, and Liu and Fagel M analyze earthquake soil-structure inter- 
actions. In the last three applications cited, the fast Fourier transform 
(FFT) algorithms were used in conjunction with the Fourier method 
to obtain numerical solutions. The development of FFT algorithms in 
the past few years has rapidly expanded the application of the Fourier 
method to include digital analysis of linear-system dynamics. 

This work was motivated by the idea that some simple and effective 
techniques based on the Fourier transform method and response- 
bound relations of linear systems can be developed for the shock and 
vibration analysis of telephone structures. 

This paper discusses the Fourier transform method applied to the 
vibration problems of linear structures and presents some simple 
relationships that make it possible to quickly evaluate and estimate 
various response parameters. In a typical engineering vibration prob- 
lem, the primary concerns are the time and frequency aspects of the 
input and output (response) and the transfer characteristics of the 
system or structure involved. The parameters or functions of interest 
are generally the waveforms of input and output, their Fourier trans- 
forms, the peak values, the frequency characteristics of these wave- 
forms, and the effects of linear filtering by structures. Of further con- 
cern are the time and frequency distributions of energy or power of 
these waveforms and, particularly, the analytical and empirical rela- 
tions among all the parameters and functions mentioned above. In 
view of the advance of computational techniques, some well-known 
and lesser-known input-transfer-output relations useful to shock and 
vibration analysis should be systematically presented. The purpose of 
this paper is to make these relations available, to discuss their implica- 
tions, and to demonstrate their applications. 

II. BACKGROUND AND THEORETICAL BASIS 

Applied to the engineering vibration problem, the basic concept of 
the Fourier method of analysis is very simple. We assume that the 
structure under consideration is linear, causal, and stable, and its 
property is completely defined by the associated transfer function 
pairs h{t) «-* H(ico), where H(iu) is the frequency response function, 
h(t) is the impulse response function, t is the time variable, « is the 
frequency variable, i is the complex unit, and the symbol <-» indicates 
a Fourier transform (FT) pair. When the structure is excited by an 
arbitrary (deterministic or random) input function, x(t) <-> X(uo), the 
response is given as y(t) «-> Y(iu). If X{iu>) and H(iu) are known, the 
response can be obtained by using the simple relationship, y(t) 
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= x(t)*h(t), where * indicates convolution, or by using Y(iw) 
= X(io))H(iu>) and taking the inverse FT of Y(iw). Other response 
quantities of interest can also be obtained from the basic FT pairs 
x(t) <-> X(ua), h(t) <-> H(iu>), and y(t) <-> Y(iu>) once their relations and 
the pertinent response parameter are found. This approach is efficient 
because no numerical integrations are required (as in many traditional 
methods of analysis), and it yields all needed solutions in both the 
time and frequency domains with simple discrete FT or FFT computer 
routines now generally available. However, the Fourier method of 
analysis is particularly attractive and powerful because of the availa- 
bility of some analytical input-transfer-output relations as described 
in this paper. These relations, in the form of simple equality or in- 
equality, can solve many practical problems. Some important pre- 
liminary background material, particularly the definition and proper- 
ties of an arbitrary shock function, f(t), will be briefly given before 
the main results. 

2.1 Shock function and related quantities 

A time function f(t),t(E [a, 6] for finite a > b ^ is referred to as 
a shock function. Without losing generality, it is assumed that f(t) 
begins at t = and has finite duration T, i.e., a = and b = T. The 
shock function f(t) could be the input x(t), the output y{t), or the 

impulse response function of the system h(t), and it could be either 
deterministic or random. Associated with a shock function f(t), the 
following quantities of interest are denned : 

Fourier Transform: 

F(iu) = f f(t)exp(-iwt)dt. (la) 

Running FT or Spectrum: 

F(t,ia>) = [ /( T ) exp (-iur)dr. (lb) 

Jo 

Energy: 

Ef{t) = £p(r)dr = ^ J~ \F(t, to.) ['do,. (lc) 



Total Energy: 



E f = r p(mt = ±- r \F(i»)\*d*. ad) 

./o ^tt y_ M 
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Instantaneous Power Spectrum: 



Pf (t,iu>) = | | F(t, io>) I s . (le) 



Power : 



p,(t) = ^) = ^L £ p/tt to)dw . ( if) 

Frequency Moment: 

M,, n = ^ P |w|»|F(iw)|da) n-0,1,2,**-. (lg) 

Time Moment: 

m fn = I* \t\»\f(t)\dt n = 0,l,2,---. (lh) 

./o 

Frequency Variation: 

Sf.n = ^ f |F<»>(tu)|d«, n- 0,1,2, ..-. (li) 

Time Van'afo'on: 

s/.n= f T \f ln) (t)\dt, n = 0, 1,2, ••-. (lj) 

Correlation Function: 

Rf(r) = /(0*/(0 = /"'/(0/(* - r)*. (Ik) 

Note that some of these quantities, e.g., M fl „ and m/,„, may not exist. 
Some simple bounds on f(t), F(u), and their derivatives, which will be 
used to derive response relations, can be readily obtained from the 
above equations. 

2.2 Some bounds ot shock function and related quantities 
2.2.1 Derivative bounds 

The following inequalities relating the derivatives of a shock function 
or of its FT with various signal parameters can be established. 

|/ (n) (0l ^ M n n = 0,1,2, ••-. (2a) 

I^WI^hs^fe] 1 - = 0,1,2,.... (2b) 

\F™(iu)\ ^ m„ n = 0, 1, 2, •••. (2c) 

|F<»>M| ^^"(^n;)* n = 0, 1,2, ••-. (2d) 
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These relations can be proved as follows. The nth derivatives f {n) (t) 
= d n f(t)/dt n and F*(tw) = d n F(iu)/do> n (n = 0, 1, 2, •••), can be 
written as 



and 



/ (B) (0 - K- I" W^W exp M)dw (3a) 

F<»>(&>) = /" (-*)"/(*) exp (-iufidL (3b) 

Jo 



Equations (2a) follows from (3a) and (lg), and (2c) follows from (3b) 
and (lh). Note that (2a) and (2c) are analogous bounds in terms of 
frequency and time moments, respectively. Equations (2b) and (2d) 
follow from (3a) and (3b), respectively, and Schwartz inequality. In 
(2b), ft is a constant and is denned as F(iu>) ^ f or |w| §= fi; i.e., in 
the practical sense, the shock function f(t) is assumed to be essentially 
band-limited. Although a shock function cannot be strictly band- 
limited, in practice most functions (except ideal impulses) possess a 
frequency value beyond which the Fourier spectra are negligibly small. 
Furthermore, because the sample rate can never be infinitely large, 
digitization of the time function always imposes a practical band-limit. 
Although the details are beyond the scope of this paper, various 
methods and criteria are available for the construction of equivalent 
band-limited signals. Equations (2b) and (2d) are also analogous 
bounds in terms of energy. 

2.2.2 Amplitude bounds 

1/(01 ^ (^V > f or all i, (4a) 

\F(io»)\ ^ {E f T)\ for all co, (4b) 

1/(01 ^%- n , for all*, n = 0,1,2,-.., (4c) 

I F(ita) | ^ T^n- , for all u, n = 0, 1, 2, • • • . (4d) 

l w l 

When n = 0, (4a) and (4b) are special cases of (2b) and (2d), respec- 
tively. From the relation ( — it) n f(t) «-> F (n) (iw), the following relation 
holds : 

1/(01 ^ l(«)-| £ f_ \F M (io>)\du>. (5) 

Relation (4c) is obtained from (li) and (5). Note that in (4c) no 
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band-limitedness requirement is needed. Similarly, the relation 
f in) (t) <-> (iw) n F(iw), eq. (lj), and 

\F(iu)\ ^ |(««)-»| C |/ (B) (0I* (6) 

J — 00 

lead to (4d). When (4c) [or (4d)] is used, better bounds are provided 
for small t (or w) when n is small, and for large t (or w) when n is large. 

2.2.3 Increment bounds 

Let 5/ - /(< + At) - f(t) and 8 F = F[i{u + Aco)] - F(io>), then 

|«/| ^ MiAt, (7a) 

\8p\ ^ miAw, (7c) 



and 



iM^O-TE 5 )]* *"•(¥)'■ (?d) 



Relations (7a) and (7c) are special cases (when n = 1) of (2a) and 
(2c), respectively. Equations (7b) and (7d) follow from (la) and the 
corresponding inverse FT relation, respectively, and from Schwartz 
inequality. Note that the band-limitedness requirement is imposed to 
derive (7b). Equations (7a) and (7c) are analogous increment bounds 
in terms of moments ; (7b) and (7d) are analogous bounds in terms of 
energy. Kak 7 proves (2d), (4b), and (7d), while Papoulis 8 - 9 proves (2b), 
(4a), and (7b). 

The equality conditions for the above relations can be easily identi- 
fied. For example, the equality holds for (2d) when f(t) = C(it) n 
for even n and a constant C, for (4b) when f(t) = C exp (iuj), 
when o} = nir/T, and n = 0, 1, 2, •••, and for (7d) when f(t) 
= 2(1 — sin Aw) 1 . Analogous conditions apply to their analogous 
counterrelations and can be established straightforwardly. 

2.2.4 Energy bounds 

Let f m = max 1/(0 1, occurring at t = t f , and F m = max \F(ita)\, 
occurring at w = uf. Then 

FIT- 1 ^E t ^ JIT for all «. (8) 

The upper bound is obvious from the definition of E { given in (Id); 
the lower bound follows from applying Schwartz inequality to the 
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relation 

F(ua) = f f(t) exp (-iut)dt, 

which results in 
|F(i«)| 2 ^ /" \f(t)\ 2 dt f \exp (-ia>t)\ 2 dt ^ E f T for all w. (9) 

III. RELATIONS FOR VARIOUS RESPONSE PARAMETERS 

The basic properties of a general shock function and its FT provide 
the basis for establishing important equality and inequality relations 
for various response parameters of linear systems. Let x(t) «-* X{iu>), 
h(t) +-» H(iu>), and y(t) «-+ Y(ia>) represent, respectively, the FT pair 
of the input, the transfer function, and the response; x(t), h(t), and 
y(t) are shock functions; therefore, the basic definitions and prop- 
erties are applicable. The functions x(t), h(t), and y(t) are assumed 
to possess time duration T x , Th, and T y and are band-limited 
within cox, uh, and coy. It is obvious that T y = T x + T y , and 
fir = min (fix, fi//)- Let x m = max |«(t)|, ft* = max |/i(OI> Vm 
= max | y(t) | , X m = max |X(tw)|, i7 m = max \H(iw)\, and 7 ro 
= max |F(tu)|, and t x , h, t y , u x , ^//, and wy be the corresponding 
times and frequencies these maxima occur. Note that t x , th, and t y or 
wx, o>H, and coy need not be equal. In addition to the well-known basic 
input-transfer-output relationships y(t) = x(t)*h(t) and Y(iu>) 
= X(iu)H(iu), the important relations that follow can be established. 

3.1 Total energy bounds 

The total energy of x(t), h(t), or y(t) is bounded by the peak ampli- 
tude of the associated time function from above and by the peak 
amplitude of the corresponding FT from below. 

X 2 m T- x ^E x ^ xlT x (10a) 

BIT* 1 <E h ^ }&T h (10b) 

YlTf 1 ^E y < ylT y . (10c) 

These results follow directly from the use of (8). 

3.2 Kinetic, potential, and dissipated energies 

For a simple oscillator with unit mass, a natural frequency a> , 
and a viscous damping X, the following relationships exist for the 
kinetic energy [_T(t) = ?/ 2 /2 <-> T(«)], the potential energy [F(0 
= u) 2 y 2 /2 <-> 7 (&>)], and the dissipated energy [D(t) = \u>lyy <-» Z)(w)], 
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of the system 



7(«) = -p [7(to)*y(to)] 
J>(«) = ^ [fW*F(to)l 

Z7T 



(11a) 
(lib) 
(He) 



where #(<) <-> F(io>) = (iu)Y(ita). The proof of these equations is 
obvious. Using the relationship Y{iu>) = H(io))X(iu>) with a simple 
FT routine, (11a) through (lie) allow various energy functions of 
time to derive directly from X{iui) and H(iu). 

3.3 Maximum response bounds 

Various upper and lower bounds of the maximum temporal response 
y m = max \y(t)\ can also be established. Let U denote lower bounds 
and Ui denote upper bounds of y m (i — 1, 2, • • •)» tn en the following 
relation holds : 



h £ h ^ y m ^ 



U\ ^ <u 6 

u 2 [ui , 

Ui 



(12) 



where 



h = max £E u (t)/ty, h = YJT y , 

t 

U 2 = {E x E h )\ U 3 = Sy, n /(ty) n , W 4 = 

Us = Sl Y Y m /w, U 6 = {Q r E x /irYH m , 



til = (Bfk/r)*, 



fir 

V Y{iu) 
-tt Y 



/\2irt y \, 
ui = (Q Y E h /Tr)*X m , 



and 



Us = ^lyX m H m /TT. 



The symbol V in the expression of Ui indicates the total variation. 

The lower bounds shall be proved first. Considering the running FT 
of y(t) according to (lb) and applying Schwartz inequality, the follow- 
ing can be written 



Y(t,iu>)\ 2 ^ f \y(r)\ 2 dT f |exp (-«W)| 2 di 



^ E y (t)t 

^ [max |v(r)|] 2 * 2 . 



(13) 
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Therefore, for all w and t, 

. , m^ /E„(t)y . \Y(t,iw)\ 
max 2/( T ) ^ ( -^ I ^ i — 4 — - (14a) 

and 

y OT ^max[^]^^, (14b) 

which gives the lower bounds l\ and U. Note that (E y /T y )* and 
[2£ y (<»)/<»]* are two special cases of h, and that h occurs at large (or 
small) t if the energy function E y {t) increases faster (or slower) than t. 
For a given structure, Y m is fixed; therefore, the bound h decreases 
with T v . Note also that, in (14a), the quantity max (iU [|^(<, ««)|/0 
should be the exact expression for h. However, the current expression 
Y m /T v is more convenient to use than the exact expression that requires 
the computation of the running spectrum over the entire t — u plane. 
Clearly, better lower bounds for the maximum response can be estab- 
lished on the knowledge of the total energy of the response than on the 
knowledge of the peak Fourier spectrum of the response. 

It is interesting to compare the lower bounds in (12) with the bound 
Liu 10 establishes for a specific case. For a simple linear oscillator with 
natural frequency co and viscous damping X, it has been shown that 
the maximum displacement response 

y m ^- |X(to)L-p= U, (15) 

where p = (1 — \ 2 )*u> . In this situation, the corresponding H(ib)) 
= (— o> 2 + 2iuu + co 2 ,) -1 and H m = (2\o> p) -1 occurring at u = uh 
= (1 — 2\ 2 )W Therefore, the bound h gives y m ^ \X(iu)\ umKtH / 
(2pa) T), and in (15) h gives y m ^ \X{iu>) | u = P /(p). Clearly, if damping 
is small, i.e., if u H ~ p, then h < h when 2\ca T > 1, U > l 3 when 
2Xo> T < 1, and h = h when 2\<a T = 1. 

The upper bounds, w x through Ms, will now be derived. The bound 
Wi follows from the use of (4a) and the assumption t y ^T y ; w 2 is obtained 
from applying Schwartz inequality to the relation y(t) = x(t)*h(t); 
Ui is obtained directly from (4c) ; u\ is a special case of u 3 when n = 1 
(Giardina 11 proves this differently) ; w 6 and u 8 follow from wi and the 
relation 

B. - i- r I F(fa) | 2 do, ^ ^ Yi ^ ?k xlHl; (16) 

Air J-a r ir 7r 

and w 6 and u 7 follow from U\ and the relation 

E y = ±- f Qy \X(i u )H(iu>)\>du,?gH 2 m E x or X 2 m E h . (17) 

^7T J-Q Y 
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Generally, better bounds can be obtained when more is known about 
the input, the transfer characteristics, and the output. For example, 
the poorest bound in (12) is obviously u s because the band-limit Sly and 
the peak values of spectra | X(iu>) | and \H(iw) | are all that are known. 
Closer bounds can be established if U Y and Y m are known, or, equiva- 
lents, E x and H m , or Eh and X m . Still closer bounds can be obtained if 
the energy functions E y or E x and E h are known. Knowledge of E V) in 
most cases, depends on the knowledge of y(t) or Y(iw); the exact 
value of y m can be readily obtained from this knowledge. Therefore, 
bounds depending on E v explicitly, such as h and u h may not help in 
practical applications. Note that bounds u h u 2 , and u 6 grow with T v ; 
therefore, they do not provide good results for large T v , i.e., when long 
excitations or small damping sustain the response for some time. The 
bound w 2 depends only on the knowledge of the input and the transfer 
function; therefore, it is convenient to use. Furthermore, it is inde- 
pendent of the output information, thus, it cannot be considered one 
of the better bounds, especially for very small or very large time t. 
Bounds Us, ut, and w 8 also depend mainly on the input and the transfer 
function, and the only output knowledge required is the band-limit 
frequency fir. Bounds u 3 and w 4 are generally better, yet they require 
the prior knowledge of t v , the time when the maximum response occurs. 
This is usually difficult to know without the waveform or FT of the 
response. Additional calculations are also needed for the following two 
bounds: for m 3 , the calculation of S v , n and the areas under the spectra 
of Y(iw) and its derivatives, and for w 4 , the absolute value of the total 
variation of the complex quantity Y(iu>) for the entire frequency axis. 
For Uz, small n provides better bounds for small t, and large n for 
large t. For w 4 , it provides better bounds for large t than for small t. 
In general, U\, u 3 , and U\ are better upper bounds, yet they are less 
convenient to use than others. 

Some interesting features about these response bounds are note- 
worthy. Drenick 12 derives the least favorable excitation for structures 
based on u 2 . An uncertainty condition T v Q,y ^ tt, familiar in band- 
limited and time-limited signal analysis can be arrived at from bounds 
h and u&. The bound w 4 can be used to estimate the settling time t f for 
the response y(t) to remain with a specified value e 10 . For example, if 
Fi(tai) is smoother, and thus has less variation than Y 2 (iu), then 
yi(t) takes longer than y 2 (t) to settle to the level e. This implies that, 
for two comparable response Fourier spectra, the smoother one corre- 
sponds to the response of a system with higher damping and the rugged 
one corresponds to the response with less damping. 
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3.4 Mean square response 

If x(t) and hence y(t) are random, the mean square response 
°v = ^C|y(0l 2 ] (where E denotes expectation) is of interest. Let 
<rl = E[\x(t)\ 2 2, and S x (iw) and S v (iw) be the power spectral density 
functions, and R x (t) and R v (t) be the autocorrelation functions of 
x(t) and y(t), respectively, then the following relations regarding o-„ 
can be found : 

a 2 v = B v (0) = f~ R h (r)R x (r)dr (18a) 

4 = f = «^ (18b) 

<jI^~ HI. (18c) 

In (18a), R v (t) = R h (t)*R x (t), R h (r) = f^h(t)h(t - r)dt = the sys- 
tem autocorrelation function, as denned in (Ik). Because the proofs of 
these equations are straightforward, they are omitted. The mean 
square response can be provided either by (18a) upon a simple time 
integration of the product of the autocorrelation functions of input 
and impulse response, or by (18b) upon the direct knowledge of the 
mean output energy or the expected square input Fourier spectrum. 
Equation (18c) provides a simple upper bound based on mean square 
value of input and the peak amplitude of the frequency response of 
the system. 

3.5 Relations regarding input and output time histories 

The following useful equations relate the time histories of input and 
output : 

.. . \dX(t,iu)/at 

= exp(t«0| ay( , | iu)/dt 



x(t) 

y(t) 



•(t)y(t) - ~ /"" Re [flMKftto)^ (20) 

v(t) = ± B^m (21) 

i- J • 
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where 

rrij = / t'h{t)dt, n = (wii/fllo), fti = {m\/m Q — m 2 /2), 

J— co 

ni = [rtii/mo — (1 + m )mitn2/2m — w 3 /6m ], • ■ •, etc., 

and the superscript c in (19) indicates the complex conjugate. 

Equation (19) expresses the waveform of input or output in terms 
of the corresponding running spectrum and can be used to recover it ; 
(20) gives the instantaneous value of the multiple of x(t) and y(t) as 
a simple frequency integral in terms of the transfer function and 
instantaneous power spectrum density of the input; (21) [or (22)] 
expresses the output (or the input) waveform in terms of the moments 
of the impulse response function of the system and the derivatives of 
the input (or the output). Liu 13 derives (19) and (20). Papoulis 14 uses 
Taylor's expansion to derive (21) and (22). Note that (21) and (22) 
are useful only if the two series converge rapidly; i.e., if the transfer 
function H(ua) only takes significant values and can be approximated 
by a polynomial of low order within the region of | w | ^ fiy. 

IV. SHOCK AND VIBRATION ANALYSIS OF SIMPLE STRUCTURAL SYSTEMS 

The relations presented above can be applied in the shock and vibra- 
tion analysis of structures. Since many structures can be treated as 
systems with a single mode of oscillation and many others can be 
solved by superposition of various modal responses in the generalized 
coordinates, the application can be well illustrated by the analysis of 
simple structures. We therefore consider a set of single-degree-of- 
freedom systems with natural frequency in the range of 0.05 to 10 Hz 
(most commercial and industrial structures have natural frequency in 
this range), and having a damping ratio of X = 2 percent and 10 
percent, respectively. These structures are subjected to a random 
transient load resembling a strong-motion earthquake accelerogram 
(Fig. 1) and four different types of idealized impulses resembling air or 
ground shocks (Fig. 2) . The various bounds of the relative displacement 
and acceleration responses of the structures are calculated according 
to (12) and compared with the actual maximum response y m . The 
results in Fig. 3 indicate that, for earthquake type of excitations, the 
lower bound h and the upper bound u 3 (with n = 0) are the better 
bounds; Ui and Mi also appear as reasonably good upper bounds 
although, for the acceleration response, they show wider fluctuations 
than u a . Bounds h and u 6 fluctuate with Y m and vary with the natural 
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5 10 15 20 25 30 

TIME IN SECONDS 

Fig. 1 — A random transient waveform. 

frequency of the structure, u ; therefore, both have the same degree of 
closeness to the real response. Bounds u 6 , u 8 , and h provide rough 
estimates of the response ; this is expected because little information is 
needed about the input, the output, or the structure. Bounds w 6 and 
w 8 are constant for the acceleration response (see Figs. 3a and 3b) 
because, for a fixed input, E x and X m are both constant, and the 
simple structures considered, H m = (2X) -1 (1 — X 2 ) -i , i.e., the maxi- 
mum spectrum of the transfer function associated with the relative 
acceleration response is a constant for a given damping, and is inde- 
pendent of co . Bounds u 6 and u%, for the displacement response, 
decrease with u> (see Figs. 3c and 3d) because the associated 
H m = (2Xp) -1 . Bounds w 2 and m 7 , for the acceleration response, grow 
with E\, which increases with a> . For the displacement response, u* 
appears also to be a reasonably good bound. 
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SINE IMPULSE 




EXPONENTIAL IMPULSE 




LINEAR IMPULSE 




SQUARE IMPULSE 



T t 



Fig. 2 — Impulse loads. 



The results for various types of impulse loads are generally similar 
and can be represented by Fig. 4. For impulsive types of excitations, 
we can make the following general observations: (i) h seems to be a 
very good bound for most cases (h tt y m when o> — > or when w is 
small), (it) the best upper bound is Uz (n = 0), which is very close to 
y m , particularly for the displacement response, (in) u 6 and u s , because 
they are constant for the acceleration response and decrease with w 
for displacement response, do not appear to be good upper bounds, 
(iv) for the acceleration response, u 2 and u 7 increase with u ] the rate of 
increase is reduced in high-damping cases, (v) for the displacement 
response, u b — > Wg and w 7 — > wi, when o> — > 0, and (vi) next to uz, u% 
and w 4 seem adequate as upper bounds. 

V. APPLICATIONS TO TELEPHONE STRUCTURES IN AN 
EARTHQUAKE AREA 

The above analysis can be applied to shock and vibration problems 
of telephone structures. Examples are the dynamic response analysis 
and aseismic design of various telephone building and equipment 
systems located in an earthquake area. Some of such applications will 
be described below. 

5.1 Aseismic design of battery stands 

It has been recognized 16 that, in large areas of the country susceptible 
to earthquakes, battery stands on floors of multistory telephone build- 
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ings require special physical design consideration to protect the 
safety of operating personnel and the reliability of telephone plant. 
Figure 5 shows a typical soft-site stand modified for service in earth- 
quake areas. It consists of bracing to support the stand laterally from 
an overhead ironwork structure and to support the cells on the shelves. 
Clearance is provided around the entire cell to allow it to rattle around 
freely in response to ground motion. The failure of the battery stand 
generally occurs from excessive acceleration, which could lead to 
sliding and rattling of cells and their banging against one another and 
against the side racks, and from excessive displacement of the stand 
that could result in overstress of the column members of the stand. 
Therefore, the displacement and acceleration responses of the stand 
are of crucial concern to its aseismic adequacy. These response parame- 
ters can be estimated very conveniently by the response-bound rela- 
tions derived in Section III. We assume that all that is known about 
the physical property of the battery stand is that it has a damping 
ratio about 2 percent and a natural frequency in the range of 3 to 5 
Hz. 16 It is reasonable to assume that the floor motion corresponding to 
a large earthquake resembles the waveform shown in Fig. 1. The peak 
acceleration amplitude would be in the range of 0.75 g because of the 
soil and building amplification of the ground motion. From the analysis 
in Section IV, the response bounds of the battery stands shown in 
Fig. 6 can be readily established. From the acceleration lower bound 
h in Fig. 6 that provides an underestimate of the actual responses, it is 
immediately clear that the natural frequency of the battery stand in 
the range of 1.5 to 3 would not be desirable. Because the high accelera- 
tion level in this range would cause lateral forces exceeding the friction 
force between the base of the cells and the supporting shelves, the 
battery would then slide and rattle violently in response to ground 
motion. The displacement bounds u 3 and h indicate a continuously 
decreasing trend with the natural frequency of the stand. The upper 
bound u 3 , which provides a conservative estimate, indicates that, if 
the battery stand is designed to have a natural frequency higher than 
5 Hz, the actual maximum displacement response would be less than 
the acceptable limit, e.g., about 3 inches. We can therefore conclude 
from such a simple analysis that, from consideration of both the 
acceleration and displacement responses, it is desirable to increase the 
frequency of the battery stand up to about 5 Hz or above. The increase 
of the frequency or, equivalently, the stiffness of the stand can be 
generally achieved by adding bracings to the wall, columns, or to 
superstructures such as the auxiliary frames. 
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Fig. 3 — Response bounds of single-degree-of-freedom structures subjected to 
random load. 
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Fig. 4- — Response bounds of single-degree-of-freedom structures subjected to sine 
impulse. 
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Fig. 4 (continued). 
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Fig. 5 — Soft-site battery stand modified for service in earthquake area. 
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Fig. 6 — Response bounds of battery stands. 
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5.2 Earthquake adequacy of community dial offices 

A community dial office (CDO) is a small telephone office housing 
mostly step-by-step equipment servicing a local community. The CDO 
buildings are generally one-story concrete block constructions, with a 
few steel-frame constructions as exceptions. 

Tests and analyses relating to typical CDO facilities 1718 have in- 
dicated seismic weaknesses in the equipment connective system, par- 
ticularly at the junction of the superstructure bracing system to the 
building. Figure 7 shows a schematic cross-sectional representation of 
a typical CDO equipment support structure. Equipment frame lineups 
are generally tied together at the 9-foot level by cross-aisle cable racks 
and angle braces (indicated as lineup braces). On the common distribu- 
tion frame (CDF) side, the frames are attached to a wood batten on 
the building wall by horizontally oriented cable racks and angle braces 
(indicated as wall braces). 

While there are several possible failure modes of CDO equipment, 
the most likely and the most critical one is the failure of the lateral 
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Fig. 7 — Typical CDO support system. 
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support system during severe ground shaking. The wood batten shown 
in Fig. 7 is attached to the building wall by either lag screws into wood 
studs, lag into lead expansion shields inserted into concrete blocks, 
or toggle bolts into concrete blocks, depending on the building con- 
struction. This batten could pull away from the wall and result in the 
collapse of the equipment. To determine whether such a failure might 
occur, it is necessary to compare the axial force of the wall bracing 
development during the earthquake with the allowable force. The wall 
force is given as F = EAb/l, where E is the modulus of elasticity, A 
and L are the cross-sectional area and length of the wall brace, respec- 
tively, and 5 is the relative displacement between the two ends of the 
wall brace. From testing or analysis we assume the transfer function 
Hj(oj) for the relative displacement response 5 of a hypothetical CDO 
system as given in Fig. 8. We wish to find the earthquake-resistant 
capacity of this system, e.g., the highest ground acceleration this 
system can safely withstand. 

A typical wall brace is made of steel with E = 30 X 10 6 pounds per 
square inch, A =0.53 square inch, and L = 5.7 inches. By testing 
we found that the maximum axial force the wall brace can take without 
pulling out of lag screws or toggle bolts is 900 pounds, which corre- 
sponds to a maximum lateral relative displacement of 5 = 0.003 inch 
between the two ends of wall brace. Based on the transfer function, 
H a (co), and the results described in Section IV, we found that a con- 
servative estimate (u 3 ) of the maximum displacement corresponding 
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Fig. 8 — Transfer function of CDO system. 
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to a 0.2 g earthquake loading is approximately 0.002 inch. Assuming 
the system remains linearly independent of the input amplitude, we 
can conclude that the hypothetical CDO system under consideration 
is capable of resisting earthquake ground acceleration up to about 
0.3 g, a motion equivalent in magnitude to the El Centro 1940 earth- 
quake. Additional braces would be needed if higher acceleration level 
is expected at the site. 

VI. CONCLUSION 

Various equality and inequality input-transfer-output relations are 
derived to facilitate the application of the Fourier transform method 
in the shock and vibration analysis of linear structural systems. In 
particular, a set of lower and upper bounds of the maximum re- 
sponse is obtained. Such bounds allow various maximum response pa- 
rameters to be easily estimated with little computation effort. The 
accuracy of the estimate depends on the degree of availability of in- 
formation about the input, the transfer function, and the response of 
the system. In general, h and u 3 appear to be reasonably good lower 
and upper estimates of the maximum response, respectively. Estimates 
of response bounds of simple systems are examined when subjected to a 
random transient excitation and impulse loadings. For systems whose 
transform functions roll off quickly at high frequencies, or for inputs 
that have fast-decaying Fourier spectra, the relations presented apply 
reasonably well. Care should be taken with short-duration impulses 
that have wide spectra. Telephone structures subjected to earthquake 
excitations, e.g., a conventional battery stand and a CDO system, are 
considered to demonstrate the practical applications. 
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